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ABSTRACT 

Stellar perturbations affect planet-formation in binary systems. Recent studies show that the planet- 
formation stage of mutual accretion of km-sized planetesimals is most sensitive to binary effects. In 
this paper, the condition for planetesimal accretion is investigated around a CenB, which is believed 
to be an ideal candidate for detection of an Earth-like planet in or near its habitable zone(0.5-0.9 AU). 
A simplified scaling method is developed to estimate the accretion timescale of the planetesimals em- 
bedded in a protoplanetary disk. Twenty-four cases with different binary inclinations(iB=0, 0.1° 1.0°, 
and 10°), gas dcnsities(0.3,l,and 3 times of the Minimum Mass of Solar Nebula, MMSN hereafter), 
and with and without gas depletion, are simulated. We find: (1) re-phasing of planetesimals orbits 
is independent of gas depletion in a CenB, and it is significantly reached at 1 — 2 AU, leading to 
accretion-favorable conditions after the first ^ 10^ yrs, (2)the planetesimal collision timescale at 1-2 
AU is estimated as: T^^j ~ (1 -I- lOOzs) x 10"^ yrs, where < is < 10°, (3)disks with gas densities of 
0.3-1.0 MMSN and inclinations of 1°-10° with respect to the binary orbit, are found to be the favorable 
conditions in which planetesimals arc likely to survive and grow up to planetary embryos, (4)even for 
the accretion-favorable conditions, accretion is significantly less efficient as compared to the single-star 
case, and the time taken by accretion of km-sized planetesimals into planetary embryos or cores would 
be at least several times of T^^^ , which is probably longer than the timescale of gas depletion in such a 
close binary system. In other words, our results suggest that formation of Earth- like planets through 
accretion of km-sized planetesimals is possible in a CenB, while formation of gaseous giant planets is 
not favorable. 

Subject headings: methods: numerical — planetary systems: formation 
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I. INTRODUCTION 

Planet formation in binary systems is an important 
issue, as a major part of solar-type stars were born 
in binary or multi-stars systems (Duquennoy & Mayor 
1991; Mathieu et al. 2000; Duchene et al. 2007) and 
approximately ^ 25% of detected exoplanets arc ex- 
pected to inhabit binary star systems (Desidera & Barbi- 
eri 2007). Most of the current detected planet-bearing bi- 
nary systems are wide S-types(Eggenberger et al. 2004), 
meaning the companion star acts as a distant satellite, 
typically orbiting the inner star-planet system over ~ 
100 AU away. Nevertheless, at least three planets— 
GJ86b(Queloz et al. 2000), 7 Cephei b (Hatzes et al. 
2003), and HD41004b (Zucker et al. 2004)-are found 
in close binary systems with stellar separation of only 
^ 20 AU. In addition, our closest neighbor— a Centauri 
AB, of which no planet has yet been detected— is also a 
close binary system with stellar orbital scmimajor axis of 
23.4 AU and eccentricity of 0.52. The binary proximity 
and eccentricity can stir strong gravitational perturba- 
tions which may pose a threat to planet formation. To 
address this issue, many researchers(Barbieri et al 2002; 
Quintana 2004; Quintana & Lissauer 2006; Quintana & 
Lissauer 2007; Quintana et al. 2007; Guedes et al. 2008) 
have performed a series of N-Body simulations and found 
Earth-mass planets can form in such close binary sys- 
tems. However, all their simulations have made an as- 
sumption that the planetary disk was initially made of 
proto-planet of at least Lunar-mass. If Earth-mass plan- 
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ets can form from Lunar-mass protoplanets, the question 
of whether Lunar-mass embryos can form and remain in 
stable orbits is critical towards understanding the like- 
lihood of Earth-like planets forming in binary star sys- 
tems. 

1.1. Planetesimal Accretion in Binary Systems 

According to the current standard planet-formation 
model, protoplanets or embryos are formed by accretion 
of 1-10 km-sized (in radii) planetesimals whose own ori- 
gin is still a big puzzle to our knowledge (Goldreich & 
Ward 1973; Weidenschilfing & Cuzzi 1993; Youdin 2008). 
The accretion from planetesimals to embryos usually fol- 
lows a runaway mode in single-star systems(Grecnberg 
et al. 1978; Wctherill & Stewart 1989; Barge & Pel- 
lat 1993; Kokubo & Ida 1996, 1998, 2000; Rafikov 2003, 
2004). However, things are entirely different for the close 
binary systems. 

Accretion in binary systems may probably not follow 
the normal runaway mode or even be suppressed be- 
cause the binary perturbations stir up the relative veloc- 
ity (Ay) among the planetesimals. Recent studies have 
found that this issue is much more complicated if the bi- 
nary perturbations are coupled with gas drag forced by 
the gas disk. Marzari & SchoU (2000) found that the 
Ay is significantly reduced by the "phasing" of plan- 
etesimals' orbits ("orbital phasing" hereafter) under the 
coupling of binary perturbations and gas drag. However, 
Thcbault et al. (2006, 2008, 2009) have pointed out that 
the "orbital phasing" is size-dependent. If assuming a 
population of planetesimals with a size distribution, such 
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as a Gaussian distribution with a medium of 5 km and 
range of 1-10 km, the average AV is still too high for ac- 
cretion in most of the region of binary systems with sep- 
arations of ^ 20 AU. Paardckoopcr ct aL (2008) investi- 
gated this issue with an evolving-disk model and found 
the conditions for planetesimal accretion were more hos- 
tile than in the non-evolving symmetrical disk. To rescue 
the planet formation in such close binary systems, Xie & 
Zhou (2008, 2009) have found two individual mechanisms 
which may reduce the effects of "size dependence of or- 
bital phasing" and favor the planetesimal accretion. One 
mechanism is "orbital re-phasing" induced by gas deple- 
tion. Taking the 7 Cephei system as an example, Xie 
& Zhou (2008) have found that "orbital phasing" is sig- 
nificantly restored for all size planetesimals after several 
timcscales of gas depletion. This "orbital re-phasing" 
reduces AT^ to very low values, favoring planetesimal 
accretion. Another mechanism is small binary inclina- 
tion. Taking a CenA as an example, Xie & Zhou (2009) 
found that the AV can be significantly reduced if a small 
inclination of such as 1° — 5° is considered between the 
binary orbit and the gas disk. However, the efficiency 
of these mechanisms may be quite diverse in different 
binary systems with different stellar mass ratios, eccen- 
tricities, inclinations, and separations. For an example, 
"orbital re-phasing" mainly depends on the timescale of 
gas depletion in 7 Cephei systems as shown in Xie & 
Zhou (2008), while it is found to be nearly independent 
of gas depletion in a CcnB as shown later in this paper. 
In addition, for highly inclined binary systems(iB > 10°), 
Marzari et al. (2009a) found gas drag has little effect on 
planetesimals because of their significant vertical excur- 
sion above the gas disk plane, and the Kozai mechanism 
strongly inhibits planetesimal accretion if is larger than 
40°. 

A crucial problem that has not been addressed in all 
previous studies is: what is the timescale for planetesi- 
mal accretion in binary systems. This timescale can be 
roughly treated as the coUisional timescale Tcoi if most 
collisions lead to accretion rather than erosion. In single 
star systems, Tcoi can be estimated as 

g 1 2a < i > m 

where n, m, Rp, a,< i > and are the number density, 
mass, radii, semi-major axis, average orbital inclination 
and surface density of planetesimals in the disk, fd and 
/ice are scaling numbers of solid surface density and solid 
enhancement beyond ice line, respectively. Here, = 1 
corresponds to MMSN, i.e., = lOgcm"^ at 1 AU, and 
fice = 4.2 accounts for solid enhancement beyond the ice 
line, which usually lies in 2-3 AU for solar type stars. 
Considering the relation, 

AV ~ < e > Vk ^2 < i > Vk 

^60(^)(^)-i/2(^)i/2 ms-i (2) 

where Vk is the local Kcplerian velocity and < e > is the 
average orbital eccentricity of planetesimals, Eq(l) can 



be rewritten as 

This estimation give a typical coUisional timescale of 10^ 
yrs roughly for planetesimals with 1-10 km sizes at 1 AU. 

In binary systems, however, neither Eq(l) nor Eq(2) is 
satisfied, and thus the coUisional timescale for binary sys- 
tems [Tj^^i) is not clear. Without any knowledge of T^^j, 
we may misunderstand planetesimal accretion in binary 
systems. For an example, if T^^j ^ 10^ yrs or longer, we 
may draw wrong conclusions on the conditions of plan- 
etesimal accretion by simulations of a binary system for 
only lO'* yrs. The conditions, such as AV , in the later 
several 10^ yrs, may become either hostile for accretion 
by orbital randomization or favorable for accretion by or- 
bital rc-phasing. Therefore, an accurate estimate of T^^j 
is very crucial to determine whether planetesimal accre- 
tion is favored or suppressed in binary systems. In this 
paper, we develop a simplified method to estimate T^^^. 

1.2. Planet Formation in a Centauri B 

a Centauri A and B are main sequence stars with spec- 
tral types of G2V and KIV, masses of 1.1 Mq and 0.93 
Mq, respectively. They are bound together as a binary 
system (or a triple star system if including the M dwarf 
Proxima Centauri. which orbits the AB pair with a semi- 
major axis over 10,000 AU according to Wertheimer & 
Laughlin 2006) with relative orbital semimajor axis of 
23.4 AU and eccentricity of 0.52(Pourbaix et al. 2002). 
As our closest neighbor in space, a Cen AB has been 
monitored for over one hundred years, and no planet has 
been identified around either a CenA or a CenB by far. 
It is concluded by Endl et al. (2001) that the mass up- 
per limit of a planetary companion is 2.5 MjupUer for 
a CenA and 3.5 MjupUer for a CenB at any orbital ra- 
dius. In other words, we still cannot rule out presence of 
lower mass planets, especially an Earth- like planet. In 
fact, recent simulations (Quintana et al. 2002; Guedes 
et al. 2008) have shown Earth-like planets could form 
around a CenB, and it is the best candidate for search- 
ing for Earth-like planets. Nevertheless, all their sim- 
ulations focus on the final stage of planet formation, 
and implicitly assume a disk of planetary embryos as 
their initial condition. Thus the remaining question is 
whether these embryos can form. Thcbault ct al. (2009) 
recently addressed this problem by analyzing the condi- 
tions for planetesimal accretion. They conclude plane- 
tary embryos formation through planetesimal accretion 
seems impossible around a CenB, unless the binary sep- 
aration was wider in its initial stages. However, their 
conclusions are limited in the absolutely coplanar case, 
where the inclination between the gas disk and binary 
stellar orbit is exactly equal to zero(zB = 0) 

In this paper, we numerically reinvestigate the condi- 
tions for planetesimal accretion around a CenB. The ma- 
jor contributions will focus on: (l)testing the efficiency 
of gas-depletion-induced "orbital re-phasing" in a CenB, 
(2)dcveloping a scaling method to estimate the timescale 
of planetesimal accretion in binary systems, and (3)ex- 
tending the study of Thebault et al (2009) by including 
the effects of binary inclinations. In section 2, we present 
the model and results of our simulations. Based on our 
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results, we discuss planet formation in a CcnB in sec- 
tion 3. Finally, in section 4, we summarize our major 
conclusions. 

2. SIMULATIONS 
2.1. Models and Setup 

Using the fourth-order Hermite integrator (Kokubo et 
al. 1998), we follow the evolution of 2 x 10"' test plan- 
etesimals around a CenB, under the coupling of the bi- 
nary's gravity and gas drag force. Planetesimals with 
random eccentricities of — 10"**, random inclinations 
of — 5 X 10^^, and random radii of 1 — 10 km, are 
initially orbiting around a CenB from 0.5 AU to an as- 
sumed outer disk-edge of 2.5 AU. This outer disk-edge is 
consistent with those values derived from Artymowicz & 
Lubow (1994), Holman & Wiegert (1999), and Pichardo 
et al. (2005). Gas drag is modeled by assuming a 3- 
dimensional non-evolving axisymmetric gas disk as that 
in our previous work(Xie & Zhou 2009). Such a gas-disk 
model neglects the reaction of gas disk to the binary per- 
turbations, which can further increase planetesimals' AV^ 
and inhibit their mutual accretion(Paardekooper et al. 
2008). Recently, Marzari et al. (2009b) found gas disk is 
much less excited by the binary perturbations if the disk 
self-gravity is considered. Therefore, our axisymmetric 
gas disk model should be a reliable simplification. 

We track all test planetesimal collisions by setting each 
body an inflated radius of 5 x 10~^ (rcoi/AU)^/'^ AU, 
where rcoi is the distance from the collision to a CenB. 
This choice gives a precision of ^1.5m s~^ in every im- 
pact velocity estimate. These velocity values are then 
interpreted in terms of accreting or eroding encounters 
by comparing them to two threshold velocities ( Viow and 
Vhigh, see details in Xie & Zhou 2009). These two thresh- 
old velocities are directly related to the speciflc critical 
energy Q* of the planetesimal. Since Q* is a broad pa- 
rameter without being well constrained, we then derived 
Viow and Vhigh from a pessimistic Q* and an optimistic 
one, respectively (Thcbault & Augereau 2007). For an 
example, if AU is greater than Vhigh ^ then the collision 
probably leads to erosion rather than accretion. 

In total, 25 simulations (see the Table 1) were 
performed with the consideration of different initial 
conditions such as: gas disk density (M^), binary's 
inclination(iB), and gas depletion. These include three 
different gas density cases: low(0.3 MMSN), medium(1.0 
MMSN) and high cases (3.0 MMSN), where MMSN de- 
notes the Minimum Mass of Solar Nebula (Hayashi 1981). 
For the binary's inclination, we study four cases (^^=0, 
0.1°, 1° and 10°) since is is probably within 10° for bi- 
nary systems closer than 30-40 AU (Hale 1994). The 
timescale of gas depiction is the same as that used in 
Xie & Zhou (2008) as Ma oc {t/rdep)'^-^ with Tdep = 10^ 
yr. The "run 00" is simulated for a singe-star system, 
from which the impact rate is normalized to in the other 
24 binary cases. The evolving timescale of "run 00" is 
10^ yrs, which is long enough to find enough impacts for 
statistics. For the remaining binary cases, the evolving 
timescale is set to 2 x 10^ yrs, since is expected to 
be much longer in inclined binary cases than T^^^ . 

2.2. Estimate of T^^^: Scaling Method 



Although the coUisional timescale, T^^^, in binary sys- 
tems cannot be derived as Eq (1) and (3) as in single 
star systems, we can estimate T^^^ if we know the scaling 
factor between it and T^^^. In this paper, subscripts "S" 
and "B" , denote the properties of single star systems and 
binary systems, respectively. 

Generally, collisional timescale is inversely propor- 
tional to the total collision cross-section of the system, 
namely, 

T^,i = Cs{NRlr\ 

T^oi = Cb{NRI)-' (4) 

where N and Rp are the number and typical radii of 
planetesimals, and Cs and Cb are the coefficients for 
single-star systems and binary systems, respectively. In 
a real system, N is too large to calculate, thus we usually 
take a system with less(such as n = 10^ - lO'') particles 
but larger radii(the inflated radius r = 10^^ — lO^"' AU in 
typical cases) for simulations. It has been shown (Brahic, 
1976) that a change in the number N of particles or in 
their radius Rp affects only the collisional rate of the 
system, which is proportional to the total collision cross- 
section. Therefore, the collisional timescales of simula- 
tions tcoi (namely the reverse of impact rate riimp) have 
a scaling relation with the collisional timescales of real 
systems, such as: 

nLp = lAfoi = l/{TgiNRl/nr') = Cs'nr\ 
nLp = ^/tfoi = l/{Tf,iNRl/nr') = C^'nr^ (5) 

Combining Eq(4) and Eq(5), the scaling factor between 
Tcoi and T^^i is written as: 

^rai — '^col * {'"^fmp / ''^imp) ) (6) 

where nf^p/rii^p is called "normalized impact rate". 
Since n^^^, nf^^ can be directly read out from our simu- 
lations (runs 0-24), we therefore use Eq(3) and Eq(6) to 
estimate T^^^. 

2.3. Results 
2.3.1. radial migration and distribution 

Figure 1 shows the planetesimal number distributions 
at three different epochs for the 24 binary cases. At 
the beginning, planetesimals orbit the central star with 
a uniform distribution(1000 per 0.1 AU, and 20000 in 
total) from 0.5 to 2.5 AU. As they are subject to the 
gas drag force, under which they migrate inward, the 
distribution then changes. The major features of these 
distribution changes, as shown in flgure 1 are: 

1) For iB < 1°, their evolutions of planetesimal radial 
distributions are very similar because their inclinations 
are too small to induce significant difference in radial mi- 
gration rate. On the other hand, for iB = 10° (10° ^ 0.2 
radian, which is comparable to the binary eccentricity 
of 0.52), the planetesimal migration rate are highly in- 
creased, resulting in a significant mass loss near 1 AU. 

2) Since denser gas leads to stronger gas drag and then 
faster inward drift, planetesimal inward migrations are 
much more significant in high-mass gas disks without 
depletion than in low-mass gas disks with depletion. 
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3)In nine cases(6 low-mass gas disk cases with is < 1°, 
and 3 medium-mass gas disk cases with depletion and 
is < 1°, namely runs 1,7,13, 4,10,16, and 5,11,17), very 
few planetesimals were lost by inward drift in the habit- 
able zone(0.5-0.9 AU) and outer region (such as 1-2 AU). 
There remains sufficient material for planet formation. 

2.3.2. normalized impact rate: n^rnp/'nfmp 

Figure 2 shows the normalized impact rate(n^p/n^p) 
as a function of radial distance to the central star (a 
CenB). The entire 2 x 10^ yrs evolution is divided into 
3 phases. As shown in figure 2, the major features are 
summarized as follows. 

Comparing different cases of different binary 
inclinations(iB = 0, 0.1°, 1.0°, 10°, from left to right, 
column by column), we find: 

1) For the coplanar case (iB=0), the impact rate in bi- 
nary systems is highly increased to as much as 1-2 orders 
of magnitude of the value in the single star case. This 
result has also been observed in figure 3 of Xic & Zhou 
(2009) and is explained as a result of the increase in the 
disk volume density induced by the progressive damping 
of planetesimals' orbital inclinations. 

2) As the binary orbit becomes more and more inclined, 
the impact rates decrease progressively. And there ap- 
pears to be a good trend: as is is increased by an order of 
magnitude, the impact rate is reduced correspondingly. 
For an example, at 1.5 AU, the normalized impact rates 
in the cases of is = 0, 0.1°, 1.0°, 10° are, respectively: 
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Focusing on the evolution of the impact rates during 
these 3 episodes(i = - 7 x lO'' yrs, t = 7 - 14 x lO'' yrs, 
and t = 14 - 20 X lO'' yrs ), we find: 

3) In cases such as runs 3, 21, 24, impact rates are sig- 
nificantly reduced with time because of the loss of most 
of their planetesimals caused by the inward drift. 

4) Despite similar mass loss by inward drift in runs 3, 
9, 15, impact rates in the later two inclined runs do not 
drop as much as in the coplanar run 3. Furthermore, 
in most inclined cases(runs 8,12, 14, 18, for examples), 
although mass losses are significant, impact rates are not 
reduced or even increased slightly. 

The reason for result 4) can be found in figure 3, which 
shows evolution of planetesimal orbital elements (i, 51) 
in cases of 14 and 17. At the beginning, planetesimals 
are excited to orbits with different i and il, resulting in 
reduction in the volume density of the planetesimal disk. 
Then, their orbits have a progressive re-phasing process 
due to gas damping. This re-phasing of i and Vl increases 
the disk's volume density, and then maintains the impact 
rate as long as the mass lose by gas drag is not too much. 
For coplanar binary case with is = 0, there is no "re- 
phasing" of planetesimals' i and 51, thus no maintenance 
or enhancement in impact rate. 



2.3.3. accretion ratio 

We define a possible accreting collision if its collision 
velocity is less than V/iig/i(see the appendix of Xie & 
Zhou, 2009 for details). Figure 4 shows the fraction 
of possible accreting collisions(hereafter, accretion ratio) 
during three episodes for the 24 runs. The major results 
are summarized as follows. 

1) During the first episode (t=0-7xl0^ yrs), a)for the 
coplanar cases(runs 1-6), the accreting collisions only 
weigh ~10-20% because planetesimals' orbital differen- 
tial phasing increases their relative velocities(Thebault 
et al. 2009); b)for the small inclined cases(runs 7-18) 
with 15=0.1° and 1°, the weights of accreting collisions 
increase to ~40-50% because small binary inclination re- 
duces the relative velocities by separating the orbits of 
bodies with different sizes(Xie & Zhou 2009); c)for the 
large inclined cases(runs 19-24) with is = 10°, the accre- 
tion weights reduce to '^30% because is=10° is sufficient 
large to pump up considerable planetesimals' orbital in- 
clinations, which give an extra contribution to increase 
the relative velocities. 

2) For the following two episodes(7 x 10"^ < t < 20 x lO'' 
yrs), we observe an evident increase in the accretion ra- 
tio between 1-2 AU in most cases independent of gas 
depletion. Most of these increases can be up to 80%, or 
even close to 100% in some cases during the last episode: 
t=14-20xl0^ yrs. However, for the two exceptions, runs 
19 and 22, we have not observed significant increase in 
accretion ratio. Their gas density is only 0.3 MMSN 
and the binary orbital inclination is as large as 10°. In 
such conditions, planetesimals are pumped up to highly 
inclined orbits with little time to experience a consider- 
able gas drag force. These cases become close to a gas- 
free case(see some similar cases in Marzari et al 2009a) 
where conditions usually become accretion-hostile be- 
cause of the progressive orbital randomization(Thcbault 
et al. 2006). 

3. DISCUSSIONS 

3.1. orbital re-phasing 

As shown in figures 3 and 5, the orbital re-phasing is 
independent of gas depiction in a CenB, and thus it is 
different from the one identified by Xie & Zhou (2008) in 
7 Cephei system. To understand the origins and differ- 
ences of these two kinds of orbital re-phasing, we should 
review the dynamics of a planetesimal embedded in the 
gas disk in binary systems (see details in Mazari & Scholl 
2000, Thebauh et al. 2006, and Xie & Zhou 2008). Gen- 
erally, the planetesimal's orbital elements follow damped 
oscillations. In figure 6, for example, the planetesimal's 
eccentricity follows an oscillation with a frequency of ^ 
(see Eq.(4) of Thebault et al. 2006) under the secular 
perturbations of the companion star. Then the timescale 
of secular oscillation, Tg, can be expressed as 

Te = 2tt/ii 

~3^AU^ ^AU^ ^Mo^ ^Mq^ 

X [l + 32(a/aB)^(l-e|)-^(^)]-i yrs. (8) 
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Because of the gas drag damping, it is also a damped 
oscillation with a damping timescale T^, which can be 
expressed as(see the appendix for details) 

Here a and as are the semi-major axes of planetesimals 
and the binary star, Ma, Mb are the masses of the pri- 
mary and companion, respectively, es is the orbital ec- 
centricity of the binary, and fg is the scaling number of 
gas disk density with respect to MMSN. Note that figure 
6 is just the solution based on secular perturbation the- 
ory (numerical solution of Eq(15) and Eq(16) in Marzari 
& SchoU 2000). Although this solution ignores the ef- 
fects of higher order and short period perturbations, it 
basically agrees with the direct N-body calculations as 
shown in figures 3 and 5. Besides, as this solution only 
considers the secular perturbations, its solution curves in 
figure 6 are much more clear than that in figures 3 and 
5, thus it is suitable for us to recognize and understand 
orbital re-phasing. 

As shown in figure 6, oscillation is over-damped in the 
inner region where gas is denser and thus stronger damp- 
ing, whereas it is under-damped in the outer region. The 
critical distance(ac) to separate these two damped oscil- 
lations can be roughly estimated by equating (ignore 
the high order terms) and T^, from which we have 

[esil - eir/'H^r/^' AU. (10) 

Although Eq(lO) cannot give an accurate value of Cc, it 
shows the dependences of Oc on as, ^b, fg, Ma, and Mb- 
Generally, under-damped oscillation is more favored in 
binary systems with high mass ratio(MB/Myi) and low 
density of gas disk. 

A reliable Oc can be derived from figure 6 by measuring 
the position where the oscillation begins to over-damp, 
namely where the equilibrium eccentricity begins to de- 
part from the forced eccentricity. The vertical dashed 
lines in figure 6 denote the location of Uc ~ 1.6 AU for 
a CenB and of Uc — 3.6 AU for 7 Cephei. Beyond the 
Oc, planetesimals with Rp > 1 km follow under-damped 
oscillation and eventually converge their eccentricities to 
e/ (forced eccentricity). This orbital re-phasing process 
only needs progressive damping, and it is independent of 
gas depletion. Within the Oc, planetesimals with i?p > 1 
begin to follow over-damped oscillation and force their 
eccentricities to different equilibriums depending on their 
radii. In this region, orbital re-phasing will be never 
reached unless the gas density depletes to some low val- 
ues. In 7 Cephei case, over-damped oscillation dominates 
the whole planetesimals disk, therefore orbital re-phasing 
require gas depletion as shown in Xie & Zhou (2008). On 
the other hand, in a CenB case, under-damped oscilla- 
tion is significant beyond 1.5 AU, therefore the orbital re- 
phasing is independent of gas depletion as seen in figures 
3 and 5. It is worthy to point out that the over-damped 
region in 7 Cephei would shrink to Oc ~ 2 AU if assum- 
ing a low-mass gas disk of only 1 MMSN rather than the 
massive gas disk of 10 MMSN adopted by Xie & Zhou 



(2008), which matches very well with the dependence of 
Oc on fg as shown in Eq(lO). 

3.2. mass left for planet formation 

Our results (figures 3 and 5) show the efiicient or- 
bital re-phasing of all planetesimals arc reached near 1.5 
AU of a CenB after ^ 10^ yrs, leading to a "possible 
planetesimal-accretion zone" there. This encouraging re- 
sult, however, may be challenged by several problems. 
The first one is, as pointed out by Thebault et al. (2009), 
most planetesimals in the "possible accretion zone" may 
be removed as a result of inward drift induced by gas 
drag if the accretion-favorable episode comes too late. 
Their results shows no object smaller than -^4 km is left 
beyond 0.8 AU after 2 x 10^ yrs. This issue, in fact, 
is directly related to the size distribution of the initial 
planetesimal population, which is poorly constrained by 
current knowledge on the planetesimal formation. In this 
paper, with the assumption of a random size distribution 
from 1 to 10 km, the inward drift problem is less severe 
or even ignorable in some cases, including runs 1, 4, 5, 
7, 10, 11, 13, 16, and 17. Furthermore, the outer edge of 
planetesimal disk we considered (2.5 AU) is larger than 
that(1.5AU) of Thebault et al.(2009). Therefore, in the 
possible accretion region, say near 1.5 AU, the mass lost 
by the inward drift of local planetesimals can be replen- 
ished by the planetesimals from outer disk(see the left 
panel of figure 7) . In addition, if the ice line is considered 
at 2 AU, the supplement of planetesimals is very consid- 
erable because of the enhancement of solid mass beyond 
the snow line (see the right panel of figure 7). In such a 
case, gas inward drift does not induce a mass loss prob- 
lem but a mass enhancement favoring the planetesimal 
accretion at 1- 2 AU of a CenB. For these considerations, 
we conclude, mass loss by the inward drift is negligible if 
Md < 1 MMSN and is < 10°. 

3.3. accretion-hostile transition 

The second problem is that even when mass is pre- 
served in the possible accretion zone, say 1-2 AU, the 
accretion-favorable conditions shown in figure 4 are only 
reached after an accretion-hostile transition duration of 
at least several 10^ yrs. The question that remains, as 
pointed out by Thebault et al. (2009), is how these 
objects might survive this long erosion period. Dur- 
ing this period, it is likely most large planetesimals will 
be fragmented into small debris which will be quickly 
removed by the inward drift because of their small 
sizes. To understand the details of this erosion pro- 
cess one needs a N-body code including fragmentation 
which is out of the scope of this paper. Nevertheless, 
we can reach some qualitative conclusions by compar- 
ing the collisional timescale(r^^;) with the duration of 
the accretion-hostile transition(T'a/j(). The argument is 
the following: if the disk is fragmented into small de- 
bris, every planetesimal should have experienced at least 
one erosive collision, or in other words, the duration of 
the accretion-hostile transition should be longer than the 
collisional timescale (Ta/it > T^i). On the other hand, if 

> Taht, then the planetesimals are able to survive 
the accretion-hostile transition . As shown in figure 4, 

should be at least 10^ yrs to survive the accretion- 
hostile transition. 



6 



According to Eq.(l), the collisional timescale in sin- 
gle star systems is about T^^i = 10"* — 10^ yrs for 1-10 
km-sized planetesimals at 1.5 AU where the center of 
possible accretion zone is located(see figure 4). If we 
take a medium value of T^^; = 5 x lO'' yrs, then adopt- 
ing the scaling factors of Eq(6) and Eq(7), the collisional 
timescales at 1.5 AU for different binary cases are roughly 
estimated as: 
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^col " 


- 10^ 


yrs. 


if 


«s = 


= 0.0, 


rpB 

^col 


- 10^ 


yrs. 
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- 10^ 


yrs. 
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rpB 


- 10^ 


yrs. 


if 
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: 10°, 



or expressed as a simple formula as: 
Tcoi ~ (1 + lOOis) X 10^ yrs, if < is < 10°. (12) 

Therefore, for the cases with ib > 1.0° (runs 13-24), 
collisional timescale is long enough for planetesimals to 
survive the accretion-hostile transition. Planetesimals 
are probably able to accrete into planetary embryos for 
the cases with 1° < < 10°, while cases with is > 10° 
are less favorable for planetesimal accretion, even though 
their Tcoi are long enough to survive the accretion- hostile 
phase. The reasons are the following: a)thc accretion 
timescale is so long(rcoi > 10^ yrs if is > 10°) that al- 
most all the planetesimals would be removed before they 
have accreted into bigger ones(see cases 20, 21, 23, 24 
in figure 1), and b) i_b > 10° would lead systems close 
to gas-free case, in which conditions become accretion- 
hostilc(see cases 19, 22 in figure 4). 

For the coplanar cases(iB = 0), as shown in figure 
4(runs 1-6), more than 80% collisions lead to erosion at 
the first several 10^ yrs. This "erosion-dominated phase" 
is longer than the collisional timescale by 1-2 orders of 
magnitude, thus planetesimal disk will be efficiently frag- 
mented into debris before conditions to become favorable 
for planetesimal accretion. This result is consistent with 
the conclusion of Thebault et al. (2009) 

For the small inclined cases(iB = 0.1°, runs 7-12), dur- 
ing the phase of the first several 10^ yrs(l — 7 x 10'' yrs as 
in figure 4) the erosion collision still weighs a large part 
(~ 60%), but not as dominant as in the coplanar cases. 
The chances for accretion and erosion are nearly equal 
during this period, and thus the final result of the colli- 
sional evolution during this phase is uncertain. This "un- 
certain phase" dominates the fate (accretion or erosion) 
of the planetesimal disk, since its duration (^ 10"" yrs) 
is much longer than the collisional timescale (Tcoz 10^ 
yrs, if is = 0.1°). For this reason, the cascs(iB=0.1°) are 
marginal cases, for which we cannot draw any conclusion 
about whether planetesimals arc possible to accrete into 
planetary embryos or not. 

3.4. planetesimals grow slowly 

As discussed above, planetesimal accretion is possi- 
ble, most probably at 1-2 AU around a CenB if 1° < 
Ib < 10°. Nevertheless, the accretion process in a CenB 
will be never like the normal accretion around a single 
star(type I runaway mode), since most relative veloci- 
ties among planetesimals are pumped up to values larger 



than their escape velocities. The accretion may proba- 
bly follow a type II runaway mode which is characterized 
with initially a slow orderly growth and then switches to 
runaway mode whenever the biggest body's escape veloc- 
ity is larger than AV^(Kortenkamp et al. 2001; Thebault 
et al. 2006, 2009; Xie & Zhou 2009). However, even 
for the type II runaway mode, planetesimal accretion, 
is significantly slowed down compared to the normal ac- 
cretion or type I runaway accretion. As shown in figure 
2 of Kortcnkamp et al. 2001, the planetesimal accre- 
tion timescales are ~ 10^ yrs and ~ 5 x 10^ yrs at 1.5 
and 2 AU, respectively. Although these results from Ko- 
rtcnkamp et al. 2001 are based on the Solar- Jupiter sys- 
tem, they have shown a similar mechanism(the type II 
runaway mode) which can be applied to a CenB. 

Besides, accretion process is further slowed down by 
binary inclination, which stirs the planetesimal disk to a 
larger scale height with lower volume density, and thus 
leading to lower impact rate. For the possible accre- 
tion favorable cases(iB = 1° — 10°) discussed in the 
last subsection, the collisional timescale is as long as 
Tcoi ^ - yrs, i-e., it takes at least 10^ - 10^ 
yrs to form planetary embryos or planetary cores. On 
the other hand, gas disks in close binary systems may 
deplete very fast (Cieza et al. 2009). For binary systems 
such as a CenB with a separation of only ~ 20 AU, the 
gas depletion timescale can be as short as a few 10^ yrs. 
Therefore, in the case of a CenB, although planetesimal 
accretion is possible, the accretion timescale would be 
comparable or longer than the gas depletion time scale. 
In such a case, even if Earth-like planets or planetary 
cores could be formed, there would be no gas left for 
them to accrete into gaseous planet. This argument may 
probably account for the fact that no gas giant has been 
found orbiting a CenB(Guedes et al. 2008). 

4. CONCLUSIONS 

We perform a series of simulations to investigate the 
conditions for planetesimal accretion from 0.5-2.5 AU 
around a CenB. In all of the simulations, we vary 
the gas-disk density(0.3-3 MMSN), binary inclination(0- 
10°), and consider both the cases with and without 
gas depletion. Three indicators— planetesimal num- 
ber(or mass) distribution, planetesimal impact rate(or 
collisional timescale), and the fraction of accreting 
collisions— are used in our discussions to determine 
whether planetesimals can accrete into planetary em- 
bryos or planets around a CenB. Our major conclusions 
arc summarized as following: 

1) Most planetesimals are able to remain in the "possi- 
ble accretion zone" (1-2 AU around a CenB) against the 
gas drag induced inward drift if the gas density is rela- 
tively small (0.3-1.0 MMSN) and the binary inclination 
is not too large(i_B < 10°). 

2) For all the cases with and without gas depletion, 
within a timescale of ~ 10^ yrs, we have observed an or- 
bital re-phasing for planetesimals at ~ 1 — 2 AU around a 
CenB(see figures 3 and 5), leading to an accretion favor- 
able condition there(see figure 4). This orbital re-phasing 
is caused by the under-damped oscillation of planetesi- 
mals' eccentricities and inclinations, which is indepen- 
dent of gas depletion and thus different from the orbital 
re-phasing identified as in 7 Cephci (Xie & Zhou 2008). 

3) We develop a simplified scaling method to estimate 
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the coUisional timescale in binary systems, T^^^. Wc find 
T^, sensitively depends on tlie binary inclination(iB) as 
Tcoi (1 + lOOis) X 10^ yrs, if < is < 10°, at 1 - 2AU 
around a CcnB. Based on the estimate of T^^, we find, 
planetesimals are probably able to survive the accretion- 
hostile transition period in the first 10^ yrs and accrete 
into larger planetary embryos if 1° < is < 10°. 

4)Planctesimal accretion into planetary embryos in a 
CcnB, though possible, take a long timescale as much as 
10^ — 10^ yrs, which is 1-2 orders of magnitude longer 
than that in normal accretion process in single-star sys- 
tems. Therefore, the formation of gaseous giant planets, 
like Jupiter and Saturn, are not favored in a CenB, since 
gas depletes very fast in close binary systems(Cieza et al. 
2009) in a timescale as short as 10^ yrs. 

In summary, although planetesimal accretion in a 
CenB is significantly less efficient and slowed-down as 
compared to single star systems, it is still possible if 
gas density is 0.3-1.0 MMSN and binary inclination is 
1° < is < 10°. These accretion favorable conditions, in 
fact, are typical values for initial gas density (Andrews & 
Williams 2005) and binary inclination(Hale 1994, Jensen 
et al. 2004, Monin et al. 2004, 2006) from current obser- 
vations. Our results support recent work by Guedcs et 
al. (2008), which has shown Earth-mass planets can be 
formed near the habitable zone(0.5-0.9 AU) of a CenB 



if the disk is initially composed of lunar-mass planetary- 
embryos. The possible accretion zone shown in this pa- 
per is roughly between 1-2 AU, which matches well with 
their planet formation zone(^0. 5-2.0 AU, as shown in 
figures 1 and 2 of Guedes et al. 2008). In addition, 
at the time of writing this paper, we note a promising 
result from Payne et al. (2009) that Earth-like planes 
can also form in the habitable zone of a CenB-like bi- 
nary systems through outward migration from the inner 
accretion-unperturbered zone(within ~ 0.7 AU). There- 
fore, by combining these studies(Guedcs et al. 2008; 
Payne et al. 2009) and our simulations, it is quite pos- 
sible that a habitable Earth-like planet may be hidden 
around a CenB. 
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APPENDIX 
GAS DAMPING TIMESCALE: To 

According to Adachi, Hayashi, & Nakazawa (1976), the average decay of planetesimal eccentricity caused by gas 
drag can be given as 

e dt 



(Al) 



where r] ^ 0.002(a/A[/)^, which describe the sub-Keplerian motion of gas, and e, i are the orbital eccentricity and 
inclination of the planetesimal. Taero is the characteristic timescale, which is related to the planetesimal size and gas 
density. With MMSN modeling the gas disk, Taero can be typically expressed as 



-If Rp 



Ma 



-1/2 



AU' 



'Mr: 



yrs. 



(A2) 



where Rp and a are the radius and semi-major axis of the planetesimal. Ma is the mass of primary star, and fg is 
the scaling number of gas density with respect to MMSN. In the binary case of our interest (slightly inclined binary 
systems), e is much larger than i and 77, thus the time scale for damping e is 



dt e 



(A3) 



As long as eccentricity is dominated by secular perturbations, we can take e as an average of forced eccentricity (see 
Thebauh et al. 2006) as 

^forced — , 9 : \^^} 

4 as I -el. 



then we finally get the damping timescale 



as .(l-e|) MA^_l 



/2 



•Mp 



yrs. 



(A5) 



Note this formula will be invalid if es is close to zero, since the evolution is then dominated by others, such as short 
period perturbations, rather than by secular perturbations. 
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Fig. 1. — Planetesimal number distribution (N/(0.1 AU)) vs radial distance (AU) to the central star. Trapezoid: t: 
t=3xl0'' yr, dashed: t=10xlO* yr, and solid: t=18xl0'' yr 



,=0, dashcd-dot: 



11 



i_B = 



c 0-100 r 
0-010 r 

0.001 L 



= 0.1 deg 

run= 1 2 



deg 



i_B = 10 deg 

run = 24 





0-100 r 
0.010 r 

0.001 L 








0.100 
0,010 
0.001 



0.100 
0.010 
0.001 








radial distance (AU) 



radial distance (AU) 



radial distance (AU) 



o 




d 



o 



o a. 



Fig. 2. — Normalized impact 
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Fig. 3. — Planetcsimal orbital elements (i, Q.) vs radial distance (AU) at 4 different epochs; 1, 6, 12, 18 xlO*^ yr for the case 14 with gas 
depletion(up 2 rows of panels) and case 17 without gas depletion (bottom 2 rows of panels). 
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dot: eccentricity of a planetesimal with iip = 1 km, and dashed-dot: = 10 km. The vertical dashed line denotes the location of ac, 
which separate over-damped and under-damped regions. 
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Fig. 7. — Number of Planetesimals within a = 1.4 — 1.6 AU vs time for run 17 (is = 1°, Md = IMMSN, including gas depletion). Left 
panel: case without considering the solid enhancement beyond ice line. Right panel: case including the solid enhancement of fg = 4.2 
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